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energy processes by taking into account the off-shellness (virtuality) 
£C) . of partons bound in the nucleon. The effect of parton off-shellness 

in deep inelastic ep — * eX scattering (DIS) and the Drell-Yan pro- 
cess (pp — > 11 X) is examined. Assuming factorization and a single- 
parameter Breit-Wigner form for the parton spectral function, we de- 
velop a model to calculate the corresponding off-shell cross sections. 
Allowing for a finite parton width ~ 100 MeV, we reproduce the data 
of both DIS and the triple differential Drell-Yan cross section with- 
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Qh. out an additional K-factor. The results are compared to those from 

perturbative QCD and the intrinsic-^ approach. 



X 1 Introduction 



One of the major goals of present day research is to study the structure 
of the nucleon and other hadrons in terms of the fundamental quark-gluon 
dynamics. In high energy hadronic processes like DIS, the Drell-Yan process, 
jet production, etc., the soft and hard subprocesses can be disentangled. The 
hard cross section can be calculated using the well established methods of 
perturbative QCD. This procedure allows one to extract the information 
about the nonperturbative quark and gluon properties in a bound state from 
the experimental data. 
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The described method, based on factorization, is analogous to the Plane 
Wave Impulse Approximation (PWIA) for the description of quasi-elastic 
(e, e'p) scattering in nuclear physics. The approximation of quasi-free con- 
stituents is valid when the binding energy is small compared to the energy 
transfered. In the theory of nuclei, the effects beyond the PWIA (such as 
photon radiation, initial state interaction (ISI) and final state interaction) 
are known to be essential for understanding semi-exclusive observables. Mea- 
surements, in which energy and momentum of the nucleon can be determined 
from the final state kinematics, offer an opportunity to study these effects 
and thus probe the nucleon interaction in nuclei PJ E] . One would like to 
gain an understanding of the hadron structure which is as good as the present 
understanding of the compositeness of the nucleus in terms of nucleons and 
their interaction. 

In the present paper, deep inelastic scattering and the Drell-Yan pair 
production are considered. Our aim is to investigate a kinematical region 
where standard perturbative QCD no longer works and where we thus need 
to model nonperturbative effects. Higher twists, suppressed by inverse powers 
of the hard scale Q 2 , are important in description of low xsj DIS 3], hadron- 
hadron collisions [3], and semi-inclusive DIS at moderate energies jS]. In 
the case of fully inclusive DIS, the factorization of higher twist contribution 
in terms of a hard coefficient and the matrix element of quark and gluon 
fields in the nucleon was proven Coefficients of the twist expansion were 
calculated in [S|. But the matrix element is a non-perturbative object and 
has to be modelled. In the present work we model the power corrections by 
dressing the active parton lines with spectral functions. Figure Q shows a 
handbag graph with the relevant initial state interactions which could build 
up a finite parton width. 

The initial and final state parton interaction effects on the observable 
hard scattering cross sections have recently attracted a lot of attention. The 
essential role of final state interactions in the interpretation of the measured 
DIS structure functions at finite Q 2 has been stressed in 1310 IHj- in several 
other calculations [311011111112111211111, non-collinear kinematics, i.e., non- 
vanishing primordial transverse momenta of the partons in the nucleon, was 
considered. The authors of [121 EE] pointed out that one gluon exchange in 
the initial state can produce a large effect in irp scattering in the framework 
of a quark-diquark model. On the other hand, parton off-shellness effects in 
DIS and the Drell-Yan process can have the same order of magnitude as those 
of the intrinsic transverse momentum [TZj. Thus, a consistent treatment of 
both is necessary to go beyond the PWIA in these reactions. In this paper, 
we develop the formalism to study these effects and apply it to calculate the 
cross sections of DIS and the Drell-Yan process. 
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Figure 1: The handbag graph for DIS and the relevant initial state interac- 
tions that could build up a finite parton width. 



Since the triple differential Drell-Yan cross section is a more exclusive 
observable than the DIS cross section, it is expected to be more sensitive 
to the ISI. The results of our calculations confirm this intuitive expecta- 
tion. By taking into account both the finite width generated by ISI and 
the non-collinearity of partons, we reproduce the experimentally measured 
fully inclusive DIS cross section and the triple-differential cross section of the 
Drell-Yan process very well. Our success in reproducing the transverse mo- 
mentum distribution of the Drell-Yan lepton pairs is particularly interesting, 
since other models disagree with experiment. In leading order of perturbative 
QCD, a delta function at zero transverse momentum is predicted. Only after 
the resummation of all orders in as, the pQCD predictions for the Drell-Yan 
pair p^-distributions are reliable ^Sl IH3 120]- Among the next-to-leading 
order contributions to the dilepton cross section, only the gluon Compton 
scattering can give non- vanishing transverse momentum (pr)- This process, 
however, contributes only in the region of very high Pt- Pt > V~M, where M 
is the mass of Drell-Yan pair. In contrast, the major part of the measured 
pairs lies in the interval < < \fM and is not described by leading twist 
perturbative QCD [21]. At the same time, none of the phenomenological 
models, including intrinsic kr approach, is able to reproduce simultaneously 
magnitude and shape of the experimentally observed distribution. In con- 
trast, as we will demonstrate below, the data can be successfully described 
by a model which allows for off-shell partons. 

The applied technique is presented in section |2] and the obtained results 
are discussed in section El followed by a short summary. We discuss the 
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interesting question of the application of factorization to hadron scattering 
in the Target Rest Frame in the Appendix. 



2 Method 

The basic tool in the calculation of hard processes is the factorization into 
hard and soft physics: 

da = Y^e 2 Ji(£,p±) ® dcr(£,p ± ), (1) 

i 

where the sum runs over all relevant parton flavors, is the charge of the 
zth type of parton in units of the proton charge e. da is the elementary cross 
section for a given process, f(£,p±) are unintegrated parton distributions 
defined as [22]: 

fitPx) = ^J d 'y (Nmyh + m\N) e^6(y + ), (2) 

where £ = p + /P + is the Nachtmann variable and p and P are momenta of 
the active parton and hadron, respectively. In [31| . the factorization in the 
form was proven in the leading power of the hard scale (photons virtuality 
in DIS and the Drell-Yan process). 

Note that the unintegrated distributions do not depend on p~ due to 
S(y + ). In other words, the parton distributions measure the correlation of 
partons with equal light cone times (y + = 0). This reflects the fact that the 
structure functions, measured in the fully inclusive DIS, depend only on P + 
in the Bjorken limit. To see this, let us consider the hadron tensor measured 
in DIS [23]: 

WW?) = 4^/ dW q - y (N\J,(y)J u (0)\N)c- (3) 
In the hadron rest frame (Mn is the nucleon mass) : 

O 2 

(q+, g_, q T ) = (-M N x Bj , , 0), (4) 

As Q 2 — > oo with XBj = Q 2 '/ ' P-q finite and fixed, g_ — > oo. As a consequence, 
the integral in (JHJ) should vanish due to the fast oscillating exponent, unless 

y + - 0. (5) 

At the same time, y_ is finite and even large ?/_ can contribute to W^ u in 
some cases. To be precise, the restriction on ?/_ is 



\y-\ < l/(M N x Bj ). (6) 
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In case of a fully inclusive process, one has the following condition due to 
causality: 

y 2 = y^y + - y 2 T > y T -> 0. (7) 

Thus, DIS in the Bjorken limit is a light-cone (y 2 — > 0) dominated process 
and the hadronic part of the DIS cross section is a function of a single variable 
P + = x Bj P + . 

The factorization formula (JIJ is valid only in the scaling limit, i.e., at the 
leading power as Q 2 — > oo. On the other hand, at moderate Q 2 considerable 
^-dependent corrections might be necessary to make predictions for semi- 
exclusive observables, e.g. Drell-Yan lepton pair production cross section and 
asymmetries. In this case, we propose the following factorization ansatz: 

da = J2 e2 i9i(^P±,P~)® d °(^P±>P~)- ( 8 ) 

i 

The difference between (pQ) and (jHJ) is precisely due to off-shellness effects 
that we aim to study. Indeed, the minus component of the free parton 
momentum is fixed by the on-shell condition: 

p 2 = p + p — p\ = ^P + P~ — p\ = (9) 

(we put the current quark mass to zero). However, since the partons are 
bound in the nucleon, © no longer holds. Thus all the 4 components of the 
parton momentum are independent and the full propagator should be used. 
In this case, the cross section is calculated using a virtuality distribution 
defined by a parton spectral function j2HI2Ej- Spectral functions of quarks 
in quark matter are, for example, calculated in 

In nuclear physics, the terms off-shellness and virtuality are often inter- 
changed. The on-shell condition for the nucleon reads P 2 = M^, where M# 
is the nucleon mass in vacuum. Thus, only 3 components of the on-shell nu- 
cleon's 4-momentum are independent. In case of an interacting nucleon, P 2 is 
no longer fixed and its distribution (spectral function) is given by the details 
of the interaction. All four components of the off-shell nucleon's momentum 
are independent. Thus, a hadron is said to be off-shell, if its momentum 
squared is different from the free hadron mass, i. e. when it is virtual. 

Partons in the nucleon are always virtual. For example, in the naive 
parton model, the parton momentum squared is p 2 = (xBjP) 2 = x^jM^, 
which is usually far from the current parton mass (= in our calculations). 
We call a parton off-shell, if the parton's momentum has four independent 
components. In this case the parton off-shellness p 2 is not fixed and should be 
integrated over. This differs from the "trivial off-shellness" of parton model, 
in which the quark is virtual, but it's off-shellness is fixed (to x^M^ ). More 
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realistically, one should include the transverse motion of partons. Then, for 
a free parton, p~ = p\/p + . In our calculations, p~ is not fixed by p + and p\. 
Instead, we integrate over all kinematically allowed p~ . 

In the following, we additionally assume that the dependence of g on p~ 
factorizes from the p_i_-dependence: 

&r = Y, e ifi&P±) ® da(Z,P±,m) <g> A(m,r). (10) 

i 

In ()10p . da(£,p±, m) is the off-shell partonic cross section and m = \pp* the 
parton's off-shellness. We choose 

m,p±) = m,p±)- (ii) 

Identifying fi(C,,p±) with the usual parton distribution functions means that 
we apply a quasiparticle picture, in which all effects involving more than one 
parton are encoded in the spectral function. The latter includes a width 
caused by parton-parton interactions (see e.g. 26\ and references therein). 

In our calculations, a Breit-Wigner parametrization for the parton spec- 
tral function A(m, T) was applied: 

A(m,T) = - 2 *\ (12) 
7r m z + 4I 

The width T of partons was considered constant. We find its value by com- 
paring the calculated cross section to the data. Note that for simplicity we 
use the same A(m) for all parton types. 

The hard part, i.e., the partonic cross section, is calculated using the 
rules of perturbative Quantum Chromodynamics (pQCD). We calculated the 
pQCD differential cross section for an electron scattering off a virtual quark 
and that for the annihilation of an off-shell quark- ant iquark pair into a pair 
of dileptons (see (J37|) below). Both off-shell partonic cross sections turned 
out to be gauge invariant due to the on-shell leptons making the modification 
of the vertex by Ward's identity unnecessary. 

The analysis of the off-shell kinematics and the obtained cross sections 
are separately given below for DIS and the Drell-Yan process. The case of 
electron scattering (section 12.1)1 is simpler and serves as an introduction to 
the calculation of the Drell-Yan pair transverse momentum distribution in 
section 12.21 
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2.1 DIS 

Ignoring the off-shellness of partons, the factorization formula (JIJ for DIS 
can be written as: 

da f . t r . _ „ , t , da 



„, E^/^/M-pi)i (13) 

1 o 

(If) =^(# + i 2 )^ + i + ?), d4) 

V d *<V on-shell * 2 * 
where are the Mandelstam variables, a = e 2 /4n, the parton quantities 

are labeled with hats, and the 5-function reflects the on-shell condition on 
the parton level: 

s + u + t = Q. (15) 

Let us consider the Bjorken limit (Q 2 — > oo with xbj -fixed, where q is 
the momentum transfer, Q 2 = —q 2 ) in the rest frame of the nucleon. In this 
limit, the partonic and hadronic invariants are simply related: 

i = t, s = £s, u = £tt. (16) 

From ()15|) and ()16|) one gets the constraint 

^-7h = -^h= XB >- (17) 

The parton model cross section of DIS is obtained: 

( da \ ^ 2 f 2na 2 (s 2 + u 2 )\ 

{mjjjo'^'^'^^-i+r)- (18) 

where 

fi(xBj) = J dp ± fi(x Bj ,p ± ), 

"LO" stands for leading order of perturbative QCD, i.e. parton model. 

For finite Q 2 , the fact that the partons are off-shell can generate large cor- 
rections to the formulas fT5|) - fl~8|) . We would like to point out the important 
analogies and differences to the on-shell case: 

• The energy-momentum conservation reads (c.f. ()15|) ) 

s + u + t = m 2 , (19) 

where m 2 = p 2 denotes the virtuality of the struck parton. 
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In case of an off-shell initial quark, we find the following relation be- 
tween the partonic and hadronic variables 

t = t = -Q 2 , s = £(s - M 2 N ) + m 2 , 

u = Q 2 - £(s - M 2 N ) = £(u + M 2 N ) + Q 2 (l - -i-), (20) 

which coincides with (II 6)1 in the Bjorken limit. We choose the z-axis 
along the incoming electron. 

The hadron light cone momentum fraction carried by the struck parton 
(£ = p + /P + ) is not equal to the Bjorken xgj, unless Q 2 — > oo. The 
relation between XBj and £ is 

Q 2 - m 2 - £(M 2 + 2p ± g ± 

Relation (J2T|) yields a nonlinear equation for xsj, because q± depends 
on XBj as follows: 



One can see that q± 2 < Q 2 and that it reaches its maximum at s ^> 
Q 2 /2x. Due to (|2T)l. the ISI in DIS can be interpreted as a smearing of 
the parton momentum fraction £ around its parton model value scgj- 
In the following three cases equation (J2TJ) simplifies: 

— Taking the Bjorken limit: 

x Bj = £■ (23) 

— Neglecting the transverse momentum of the struck parton inside 
the nucleon, but keeping m 2 ^ 0: 

Q 2 

XBj = ^Q 2 -™ 2 -ttM 2 N-f hS^' 

— Taking into account both the parton's transverse momentum and 
off-shellness, but considering the limit s 3> Q 2 /2x, s 3> M^: 

Q 2 

XBj = * Q 2 - m 2 + 2\p ± \y/Q* cos(0)' (25) 

where is the azimuthal angle of the quark momentum. As Q 2 
goes to infinity, equation (J25|) coincides with (J2HJ). 
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• The off-shell partonic cross section is 

(**) = 2 -^(s 2 + u 2 )S(s + u + i-m 2 ), (26) 

\dtduj off-shell t s 

where u and s depend on m 2 via (j2D|) and (|2*T|) . 

Therefore, the leading order expression for the Lorentz invariant DIS cross 
section (JT5J) is modified by the ISI as follows: 



(27) 

off-shell 



To compare to the experiment or to the leading order cross section (|18p. 
we also need to change variables from partonic u to hadronic u or (^bj is 
related to the hadronic Mandelstam variables (s, t, u) by lfI7jl ). We choose 
the following independent variables for the hadronic cross section: 

s, t, x B j. (28) 

The partonic cross section depends on: 

s, t, u, m 2 , f , p ± . (29) 

We have related the partonic s to s by fl2()(l . The transformation from one 
set of variables to the other is done in the following way: 

( Jf ) =\ dil \-lTjA S(x Bj - x Bj (s,t,u)), (30) 
\ dtdx Bj j jgj J \ dtdu ) jgj 

where (da /dtdu) is given by (J27|) and function of the variables 

(j2Hl) is defined by (|2*Tj) . We note in passing that (da/ dtdx Bj) is negative, 
while (da /dtdu) is positive. This has to be taken into account in (}3T)j) by 
an appropriate choice of integration boundaries. From equations P^j). (|27ji. 
(|3*0)l. we obtain: 

/;2 i „-.2\ 

5 (s+u+t-m 2 ^ 5\x Bj - x B j(s,t,u,C,,m 2 ,p_ L )^ ,(31) 



o 

[s 2 + u 2 

52 



8" 
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where s = £(s — Mjy) + m 2 . The integration over u can be done using one of 
the 5-functions. The result is: 



da \ ^ 2-nae 2 



(^) Br E^AKr)/^/,((iJ 

X ( 1+ (S-M^Wp ) *(*«-*«(M,4,£,mVj), ( 32 ) 

where XBj(s,t,u,£,,m 2 ,p±) is given by (|2T|) and {t = — £ — £(s — Mjy). The 
5-function can be used to perform the integration over the azimuthal angle of 
the parton momentum. The remaining three integrations must be performed 
numerically. The limit s ^> Mfj,Q 2 fix was taken for simplicity. For the 
unintegrated parton distributions f(£,p±) we use the factorized form (}4*2)) 
discussed in more detail in the next section. The results for DIS are presented 
in Sec. 13.11 

2.2 Drell-Yan process 

We applied the same technique to calculate the cross section of the Drell- 
Yan process (pp — > X + In this case, an off-shell quark-antiquark pair 
annihilates into a pair of leptons. The virtuality of the quark (antiquark) 
coming from the target proton {m\ = p 2 ) and that of the antiquark (quark) 
coming from the projectile proton (m| = p\) are in general not equal. We 
assume, however, that their distributions A(m) are the same. 

The connection between the observables and partonic variables in case of 
two off-shell particles is more complicated. Moreover, the choice of proper 
partonic variables is frame dependent. We obtain the following kinematic 
equations in the hadron center of mass system: 



999 i (m? + p? , ) (mi + pi , ) 

M 2 =m 2 1 + m 2 2 + £i£ 2 Pf P 2 + + r/ p-p+ - 2 ^ ' ^ 



XF - ~s^W 2 - 6Pl + —JJf J ' { ) 

Here, we have used: 

tLx=Pi/PT, &=pt/P 2 + , (34) 

Pi (P2) is the 4-momentum of the target (projectile) hadron, pi^ denote 
momenta of the annihilating quark and antiquark, M 2 is the invariant mass 
squared of the produced leptons. S denotes the hadron center of mass energy 
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squared. The Feynman variable is defined as xp = p z /(pz) max , where p is 
the lepton pair momentum. In some works, an approximate definition for 
the Feynman variable is used: xf ~ 2p z /\fS. We used the exact definition 
that can be written in the hadron center of mass system as follows 



Xp = 



Pz 



(Pz 



2p z VS 
S-M 2 ' 



(35) 



Experimentally observed Drell-Yan pairs have small M 2 compared to S, so 
the difference between the two definitions of xp are small. However, in the 
Drell-Yan scaling limit, M 2 ~ 5" and the difference is finite (see formula (p)3*j) 
in Appendix lA"j) . 

One sees that the definition of £2 is analogous to the DIS case, whereas 
the target's momentum fraction is defined as a ratio of minus- components. In 
some articles, alternative definitions are used, for example £1 = p% /P\ ,£,2 = 
P2/P2 ■ The choice of the definitions is based on the behavior of the 
hadron momenta in the Drell-Yan scaling limit (S — > 00). The argument 
is presented in Appendix [X] One might prefer to do the calculations in the 
target rest frame, since the connection between the observables (S, M 2 , xf, 
Pt) on the one hand and the partonic variables m 2 , pi±) on the other is 
simpler in this case (cf. (JSTj) in Appendix EJ). However, factorization in the 
form (P"U|) is not applicable in this frame of reference (for a detailed discussion 
see Appendix |A"|) . 

We have calculated the pQCD cross section of the off-shell quark-antiquark 
annihilation into a pair of dileptons: 



da 



dp [dp ' 2 



167ie[e'2S 2 N c J (pi ■ p 2 ) 2 — m 2 m\ 



S{pi +p 2 



P' 2 ),(36) 



where p [ 2 are the three-momenta of the leptons, e [ 2 their energies, and e q 
the parton charge in units of the proton charge, color factor = 3. 

The off-shell partonic cross section, differential over the Drell-Yan process 
observables - mass M, Feynman variable Xp, and transverse momentum p-p 
of the lepton pair - is: 



da 



dp[ dp'2 

dM 2 dx F dp 2 T ~J 2e[ 2e' 2 ' 



t + u 2 — m\ — m\ + s{rri\ + m 2 ) 2 



x 5(pi + P2 - p'i - p' 2 )S{p - p[ - Pa); ( 37 ) 

(38) 



a 2 e 2 q (S - M 2 



VSEM±8Nj( Pl -p2) 2 - 



1 1 

mfm 2 
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After performing analytically the seven integrations over non-measured quan- 
tities, four (^-functions are integrated out and the remaining four preserve the 
correct relation between the hadronic and partonic variables (cf. (|3*3*j) ): 

da 

dM 2 dxFd-p\ 

M 2 / _ r, , , n\ E^ / . n n\ E 



K 



— (M 2 + (mi + m 2 ) 2 ) + — (4e 2 - m\) - - e 1 (m 2 + m 2 - 



m 2 



(m 2 - m? - ™ 2 - ^ 2 Pf p 2 + - K 7^r m i/ Pl±) + ^) 

X H f "^m^1 6P2 J 

x5((pi±+p 2 ±) 2 -^)5(E- ei -e 2 ); (39) 



M 4 iV cV /( Pl .p 2 ) 

In dm, 



«' = . q . (40) 



2 — mfm 2 . 



Cl s 1 ( fl P f + K±|L>) , £2 = i ( 6 tf + (m ^ J 



Using the ansatz for the case of two off-shell partons in the initial 
state, we obtain the hadronic cross section by integrating over the masses 
and transverse momenta of quark and antiquark: 

dcr r r r°° r 00 

dMH^d& ^r^l^k dmi l dm2 Jo d H rf & A (^)A(m 2 ) 

x/,(^6,fc)/ i (Q 2 ,& P - 2 x) 5S3 f^.(4l) 

The integration in (j¥T)) is 8- fold, da is given by equations (J3~9"|) and (jjOj) . 

The common parametrization for the unintegrated parton distributions is 

f(Q 2 ,P±,0 = f(p±)-q(Q 2 ,0, (42) 

where 

m) = j^eM-§r 2 }, (43) 

and q{Q 2 ,C) is the conventional parton distribution. For the latter, we have 
used the latest parametrization by Gliick, Reya, Vogt jHOJ - The mean pri- 
mordial transverse momentum of partons is 

< p 2 ± >= AD 2 . (44) 
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The Gaussian form of f(p±) allows the analytical evaluation of the inte- 
grals over pu_ and p2±. Then, we are left with a four- dimensional integral to 
be done numerically. In the special case of a constant spectral function width, 
one of the integrals over the off-shellness can be reduced to a superposition 
of special functions (incomplete elliptic integrals). 

We compare the result of our model, in which the partons in the proton 
have a finite width, with the experimental data and with the cross sections 
obtained by two other methods (LO pQCD and the intrinsic-fey approach). 

In fey-factorization, the formula 

do = /(Ci,Pj_i)/(6>P-L2) ® d£(fi,6,£]_i,p.L2) (45) 

is used, where da is the Born cross section for the qq annihilation into a 
pair of leptons, f(£,p±) is the unintegrated parton distribution defined in 
(J2J). A proof for the fey-factorization in the Drell-Yan process is given in the 
leading twist in pH] I32j. In this case, the primordial transverse momenta of 
the q and q have (in general, non-zero) values defined by these distributions 
in the same way as the usual integrated parton distributions define the light 
cone fractions of the parton momenta (p + for the projectile parton and p~ 
for the target parton). In j^U E2], the fourth component of the parton 
momentum (p~ for the projectile parton or p + for the target parton) is set 
to zero due to the following reason. For large hard scales M, the projectile 
parton momentum is P2 = (pjj j P2 1 P2±) ~ \ 2 , X), where A = m 2 /M. 

The parameter A is small for M > 1 GeV, since the parton off-shellness and 
transverse momentum are related to the inverse of the confinement radius and 
do not scale with M. There exist several parametrizations of unintegrated 
parton distributions f(£,p±)- 

A phenomenological 'intrinsic-fey approach' has been developed on the 
basis of the fey-factorization theorem. In this model, the unintegrated distri- 
butions are taken in the form (|4Hj) . ()44jl . An additional difference from [3*T|l3*2] 
is that the smaller light cone component of the parton momentum is put to 
its on-shell value: p 2 = 02±/pti which is small, but not zero. This approach 
is well described in the literature [HI [121 HH1 Ej and proves to be very useful 
for the calculation of cross sections and asymmetries of different processes. 
It is obtained from (f4*Tf by putting all parton masses to and dropping the 
mass integrations and spectral functions. 

In the works El the Drell-Yan process was studied in the scope 
of the intrinsic-fey approach. We can obtain the partonic Drell-Yan cross 
section in the Ay-factorization approach (with on-shell partons) by putting 
m\ = = in (jSHJ), (|4()jl . In particular, the following kinematic relations 
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Figure 2: Relative deviation of the calculated DIS cross section from eq. fl~8|) 
for the range of parton spectral function widths 5 MeV to 0.9 GeV. Mean 
intrinsic transverse momentum is 0.3 GeV. Q 2 = 1 GeV 2 , s S> Q 2 . 



are ensured by the 5-f unctions in this case: 

M 2 = fcfcPfi? + ~ 2 ^ ■ 



The authors of fHJ IS] used the parton model partonic cross section and 
the approximate kinematical relations (|H3j) for simplicity. In section \3.2\ 
we compare the cross section calculated in our model with off-shell partons 
to the result of the fc^-factorization approach. In order to perform such a 
comparison, we have calculated the Drell-Yan cross section in the intrinsic- 
kx approach by using the full on-shell partonic cross section and the exact 
kinematics (|4fij) . 
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Figure 3: Relative deviation of the calculated DIS cross section from the 
result of the on-shell calculations taking into account only intrinsic transverse 
momentum effects. Q 2 = 1 GeV 2 , s > Q 2 . 




Figure 4: Relative deviation of the calculated DIS cross section from eq. (fTHJ) 
for different Q 2 . The spectral function width is fixed to T = 0.2 GeV and the 
transverse momentum dispersion to D = 0.1 GeV, s ^>Q 2 . 
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3 Results 



3.1 DIS 

The results of our calculations for the deep inelastic electron-proton scat- 
tering cross section for a range of widths as compared to the parton model 
(eq. ()18jl ) are shown in fig. El We have found that there is a moderate ef- 
fect of the initial state interaction in DIS in the region of small Bjorken xgj 
and low momentum transfer Q 2 . The cross section deviation reaches 50% 
at Q 2 = 1 GeV 2 , if the parton spectral function width and mean transverse 
momentum are both equal to 300 MeV. In fig. one can also see that the 
cross section calculated in our model differs from the LO even when the par- 
ton width is negligibly small (5 MeV). This effect is due to the non- vanishing 
intrinsic transverse momentum. 

To separate the effects of the parton off-shellness from those of the intrin- 
sic transverse momentum, we plot the relative difference between the result 
of our model with off-shell partons and the calculations taking into account 
only the intrinsic transverse motion (fig. EJ). To obtain the cross section in 
the latter approach, we put T to zero in the formulas of section 12.11 thus 
forcing the parton on-shell. It is seen that this difference amounts to at most 
10 % of the cross section. 

For values of Q 2 above 25 GeV 2 , the initial state interaction in DIS gives 
at most a 5% deviation from the lowest order cross section (fTH)) . The Q 2 - 
suppression of parton virtuality and intrinsic transverse momentum effects in 
DIS is illustrated in fig.E] For most of the experimentally investigated values 
of Q 2 , the ambiguity in the parton distribution function parameterizations 
due to the renormalization scale uncertainty is of the same order as the ISI 
effect in DIS. 

The difference between the off-shell result and the leading order cross 
section at Q 2 > 2 GeV 2 is too small to be resolved by the present experiments. 
In the region of Q 2 < 2 GeV 2 , the difference is 30 - 40 %, which should be 
observable. However, in order to make a quantitative comparison to the 
experiment at such low Q 2 and XBj, we would have to incorporate into our 
model other effects, such as resonance production and diffractive scattering 
[HHl EU [3*5] . We conclude that, using the model described in the present 
paper, we cannot extract the value of the parton width in the nucleon from 
the DIS data. This is the result expected by the analogy to nuclear physics, 
because the DIS cross section is too inclusive. On the other hand, the DIS 
data do not contradict the assumption of a finite parton width in the proton. 
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3.2 Drell-Yan process 

In contrast to the DIS case, the effect of parton off-shellness on the transverse 
momentum distribution of the Drell-Yan lepton pairs is substantial. In this 
section, we present the Drell-Yan triple differential cross section calculated 
by the method described in Section 12^21 We compare the result of our model, 
in which the partons in the proton have a finite width, with the experimental 
data and with the cross sections obtained in two other approaches (LO pQCD 
and standard /c^-factorization). 

Calculations using LO pQCD and collinear factorization analogous to 

^ = /(6)/(6)®t»(6,6) (47) 

give a simple result for the triple differential Drell-Yan cross section (jp?- 
distribution of the dileptons) - a 5-function at zero Pt- This follows from 
the fact that the annihilating q and q in this approach are collinear with the 
corresponding hadrons, thus the qq pair has no transverse momentum in the 
hadron center of mass system. Therefore, the resulting lepton pair cannot 
gain any transverse momentum in this model. In contrast, the experimentally 
measured transverse momentum distribution of the dileptons extends to pt = 
4 GeV at a hard scale (the mass of the lepton pair M) as high as 8.7 GeV. 
Note that NLO corrections do not cure the disagreement with the data. The 
Drell-Yan pair p^-distribution obtained in fixed order pQCD is divergent at 
Pt = 0. A resummation of an infinite series of diagrams is necessary to obtain 
a finite value for the triple differential Drell-Yan cross section at pt = in 
pQCD with on-shell partons ^H]- The resummed pQCD cross section is in 
qualitative agreement with experiment [T5] . 

In order to analyze the effect of a finite parton width and distinguish it 
from the effect of the intrinsic transverse momentum, we have performed the 
calculations in both the intrinsic-fcr approach and in our model allowing for 
off-shell partons. We used the formalism developed in Sec. 12.21 to calculate 
the cross section of the Drell-Yan process in the kinematics of the experi- 
ment E866 (HEl E] in the intrinsic-fey approach and for a parton off-shellness 
distributed according to the Breit-Wigner spectral function (JT2J). 

We present the obtained cross sections for different values of the parame- 
ters in figures El and El We illustrate in fig. |3]that the slope of the distribution 
mainly depends on the dispersion of the intrinsic transverse momentum (£)), 
which is proportional to the primordial transverse momentum of the parton 
(see f!44j) ) . In the limit, in which the dispersion of the intrinsic transverse 
momentum (D) goes to zero, the leading order result of perturbative QCD, 
i.e. a sharp peak at px = 0, is restored. 
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Figure 5: Calculated distribution of the Drell-Yan lepton pair's transverse 
momentum in the /c^-factorization approach for different values of the parton 
primordial transverse momentum dispersion. M = 4 GeV, xf = 0.1, T = 0. 




Figure 6: Calculated distribution of the Drell-Yan lepton pair's transverse 
momentum in our model for different values of the parton spectral function 
width. M = 4 GeV, xp = 0.1, D = 0.7. The solid line gives the result of a 
calculation with T = 0. 
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On the other hand, the parton width variation leads to changes of the 
cross section magnitude and influences the behavior of the distribution in 
the region of low pp (see fig. IBJ). One can also see in figure EH that our model 
approaches the result of the standard intrinsic-fey method as the parton width 
(r) goes to zero. At finite width, the shape of the cross section obtained in 
our model is different from the result of the intrinsic-fey approach in the low 
Pt region. Also, the magnitude of the cross section is different. This indicates 
that some additional nonperturbative effects are included via a finite parton 
width. 

In figures I71TH21 calculations both in the model with off-shell partons and 
in the standard fcy-factorization approach are compared to the data of the 
Fermilab experiment E866 for the continuum dimuon production in pp col- 
lisions at 800 GeV incident energy. In this experiment, both the double 
differential Drell-Yan cross section da/dM 2 dxF (data published in [HE]) and 
the triple differential cross section da /dp (data published in 21) were mea- 
sured in a wide range of M and xp (p is the lepton pair's momentum). The 
Pr-distribution was obtained in terms of the triple differential cross section 
averaged over the azimuthal angle of the lepton pair 

da _ 2 da 2 r da <^j\/f 2 (43) 

dp ny/s dxpdp\ ny/s J dxpdp^dM 2 

bin 

The data points were averaged in several bins in M and xp. The xp binning 
is responsible for the wiggly structures both in the data and some of our 
calculations. 

The result of the standard /cy-factorization approach is shown in figures U\ 
and|Hl (solid line). The slope of the cross section can be reproduced by an 
appropriate choice of the single parameter (D) of the intrinsic transverse 
momentum distribution, given by fi3|) . The optimal values for D are 500 — 
600 MeV, which correspond to 

< v\ > 1//2 = 1-0 — 1.2 GeV. (49) 

A slightly smaller value for this parameter was obtained in [T3*| IT4*] from the 
analysis of the data of the experiment E744 on Drell-Yan cross section in 
pp collisions at 400 GeV incident energy: 

< P± > 1/2 = 0.8 - 1.0 GeV. (50) 

Still, the data are overestimated by a factor of 2 — 3, depending on the 
mass of the Drell-Yan pair (M). Dashed lines in figures U\ and El illustrate 
that the data can be fitted by introducing an additional overall factor (K). 
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(a) 4.2<M<5.2 GeV 




(b) 5.2<M<6.2 GeV 



Figure 7: Drell-Yan cross section in /^-factorization approach compared to 
the data for continuum dimuon production in pp collision from [21]. Varying 
the dispersion of intrinsic transverse momentum (D), the slope of the distri- 
bution is fitted (solid line). An additional overall K-factor is necessary to 
reproduce the cross section amplitude (dashed line), —0.05 <xf < 0.15. 
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Figure 9: The Drell-Yan cross section as calculated in our model (solid line) 
compared to the data of the Fermilab experiment E866 for the continuum 
dimuon production in 800 GeV proton, 4.2 < M < 5.2 GeV, -0.05 < x F < 0. 15. 
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Figure 10: Same as fig. El only for a higher mass bin: 5.2<M<6.2 GeV. 
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Figure 11: The Drell-Yan cross section as calculated in our model compared 
to the data of E866 on pp^fi+^X, 6.2<M<7.2 GeV, -0.05 < x F < 0.15, 
and to the data of E772 on pd ^ fi + n~ X , 6<M<7 GeV, 0<x F <0.3. See 
main text for more details about the different lines. 



> 

CD 

CD 

5 0.1 

Q. 

"a. 

LLJ 0.01 



1E-3- 



-I — i — | — i — | — i — | — i — | — i — | — i — | — i — | — i — p 




D=0.55 GeV, r=225 MeV 

■ Experiment E866 
♦ Experiment E772 



— i — 1 — i — 1 — i — 1 — i — 1 — i — 1 — i — 1 — i — 1 — i — 1 — i — 
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 



P T , GeV 



Figure 12: Same as fig. ^2 only for a higher mass bin: 7.2<M<8.7 GeV. 
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The discrepancy between the calculations and the data is larger for higher M. 
Thus, in the fc^-factorization approach, the magnitude of the cross section 
cannot be correctly reproduced. An additional overall .fT-factor is necessary 
that reflects the importance of higher order corrections to the Drell-Yan cross 
section. 

In contrast, the calculations with a finite parton width yield not only the 
experimentally measured shape of the cross section but also its amplitude 
without any X-factor. Instead, we take care of the higher-twist and NLO ef- 
fects assumed to be contained in K by introducing the physically transparent 
off-shellness (T). Note that we work in leading order of perturbative QCD 
concerning the processes that enter the calculation of the Drell-Yan cross sec- 
tion. At NLO of perturbative QCD, additional processes contribute, namely, 
gluon bremsstrahlung (qq — > l + l~g), gluon Compton scattering (gq — > l + l~q 
and gq — > l + l~q), and vertex corrections. The bremsstrahlung, along with 
diagrams of even higher order and with gluon exchanges between the active 
(anti-) quark and the spectators, contributes to the quark width. Therefore, 
going to higher orders of standard perturbation theory, while dressing quark 
lines with spectral functions, would be a double counting. In other words, 
we have already included part of the NLO processes by our finite width cal- 
culation. Some processes beyond the NLO are also included. On the other 
hand, our model can be improved by taking into account also gluon Comp- 
ton scattering and vertex corrections. In this case, the gluon line should be 
dressed, too, and the number of model parameters increases. In the present 
paper, we concentrate on the effect of a finite quark width on observables 
and consider only the leading reaction mechanisms, i.e. qq annihilation for 
the Drell-Yan process. 

The comparison of our results with the data is given in figures IUHT21 The 
values for the average parton primordial transverse momentum 

< pi >V2 = o.9 _ n GeV (51) 

are compatible with those existing in the literature (joHj) . Allowing for off- 
shell partons, we eliminate the need for any K-factor. Choosing V in the 
order of 100 MeV (cf. table [T] for details), both the amplitude and the slope 
of the cross section are well reproduced. 

The dependence of the optimal values for the parameters (dispersion D 
and width T) on the mass of the Drell-Yan pair was obtained by fitting 
experimental data within different bins of M independently. The result is 
presented in figures I71H21 and in table HJ Note that the varying quality of 
the data in different mass bins leads to large uncertainties in the extraction 
of the width. In table [T] we present the average values and uncertainties for 
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M 


4.2-5.2 


5.2-6.2 


6.2-7.2 


7.2-8.7 


D 


450 ± 100 


530 ± 70 


540* 


550 ± 60 


r 


65 ±20 


200 ± 75 


210* 


225 ± 75 



Table 1: Optimal parameters for different masses of the Drell-Yan pair, 
— 0.05 <xp< 0.15. All values are in MeV. Values denoted with stars are 
trend-average and not best-fit. See main text for details. 

D and V. The latter have been obtained by analyzing the \ 2 values as a 
function of D and V. We find that the optimal V increases with the hard 
scale (the mass of the Drell-Yan pair). The dependence of V on M indicates 
that, at higher scales, partons with broader spectral functions are probed. 
We did not study the dependence of our parameters onxp. 

The analysis of the data in the mass bin 6.2 < M < 7.2 GeV calls for 
more discussion. As shown in figure 11, the best fit (dashed line) to this 
data set leads to values for both parameters (D=0.43 GeV, T=l.l GeV), 
which are not in the trend set by the fits to the other three data bins (cf. 
table Q). Thus, we did not trust this fit and sought for more experimental 
input. For pt > 1 GeV, the data of the experiments E866 on pp — > [i + fi~X 
and E772 [37] on pd — > fi + {i~X agree very well in all the mass bins, except 
the one of figure ^JJ (see figure EI for example). Therefore, we compared 
our fit (dashed line) to the experimental data on the pd cross section from 
E772 in approximately the same mass range (figure ITT]) . One can see that 
the calculations with D = 0.43 GeV and T = 1.1 GeV (dashed line) do 
not reproduce the high-p^ part of the pd cross section. On the other hand, 
if the trend-average values from table [JJ are applied (D = 0.54 GeV, T = 
210 MeV, solid line in figure ITT]) . the cross section calculated in our model 
both describes the E866 data on the border of experimental error bars and 
reproduces the pd cross section of E772 at pt > 1 GeV. 

Allowing for a finite parton width and using a single-parameter form 
for the parton spectral function, we account for non-perturbative effects, 
including the K-factor. The result of the collinear factorization and fixed 
order pQCD (5-peak at px = 0) is not reached in the experiment even at 
masses of lepton pairs as high as M ~ 16 GeV. There is one area of hard 
scales, where the intrinsic kp approach seems to reproduce the cross section 
with good accuracy: at low M the X-factor of the intrinsic-fcr approach is 
closer to 1. As the dilepton mass goes higher, the measured distribution is 
getting more sharply peaked. This suggests that the result of LO pQCD 
might be recovered at Drell-Yan pair masses, which are higher than those 
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yet observed. On the other hand, our model allowing for off-shell partons 
with finite width works well for all hard scales M. 

4 Summary and outlook 

The research presented here reveals the importance of the parton initial state 
interaction for the analysis of high energy processes. We have developed a 
formalism to study the quark structure of hadrons going further than the 
widely studied picture of collinear noninteracting partons. The parton off- 
shellness effects missed in the standard treatment were taken into account by 
dressing the parton lines with spectral functions and using the factorization 
assumption. In this way higher twist corrections to standard pQCD were 
modelled. 

We have calculated the cross sections of deep inelastic ep scattering and 
the Drell-Yan process pp — > l + l~X in the model allowing for a finite par- 
ton width. Off-shellness effects arise from the fact that the time-like light- 
cone momentum of the parton (p~) is not fixed by an on-shell condition 
(p~ = p\/p + ) or by a collinearity condition (p~ = 0). Since the partons 
in the proton interact, p~ is in fact distributed with some finite width. To 
disentangle the off-shellness effects from the effect of the parton primordial 
transverse momentum, we have additionally calculated the Drell-Yan cross 
section in the standard intrinsic-/^ approach. The obtained cross sections in 
both models were compared to the data on the triple differential cross section 
of the process pp — > l + l~X. 

We have found a moderate effect of the initial state interaction in DIS 
in the region of small Bjorken Xbj and low momentum transfer Q 2 . For a 
parton width of 300 MeV, the cross section increases due to the finite quark 
width in the proton reaches 10% at Q 2 = 1 GeV 2 . On the other hand, the 
effect is Q 2 -suppressed. For values of Q 2 above 10 GeV, the initial parton 
off-shellness generates only at most 2% of the cross section. For most of 
the experimentally investigated values of Q 2 , the difference between the off- 
shell result and the leading order cross section is too small to be resolved 
by present experiments. We conclude that the value of the parton width in 
the nucleon cannot be extracted from the DIS data, because the DIS cross 
section is too inclusive. This is the result expected by the analogy to nuclear 
physics. On the other hand, the DIS data do not contradict the assumption 
of the finite parton width in the proton. 

In contrast, we discover a substantial contribution of the parton off- 
shellness to the transverse momentum distribution of the high-mass virtual 
photons produced in high-energy hadron-hadron collisions in the whole range 



26 



of hard scales, at which the cross section has been measured. The triple dif- 
ferential Drell-Yan cross section is a more exclusive observable than the DIS 
cross section. That is why the effect of the parton off-shellness was expected 
to be larger in the Drell-Yan case. Our results confirm this expectation. 

Although the intrinsic-/^ approach alone can reproduce the slope of the 
experimentally measured distribution of dileptons, an overall K-factor is nec- 
essary to fit the data. Shape and magnitude of the cross section are much 
better reproduced by a model that allows for off-shell partons. In particular, 
one can fit the data without a K factor. The parton width in the proton was 
estimated from the comparison to the data. For a mass of the Drell-Yan pair 
of 4.2 — 8.7 GeV, the best fits were obtained with quark (antiquark) width 
of 50 — 250 MeV and intrinsic transverse partonic momentum dispersion of 
400 — 600 MeV. This corresponds to a mean primordial transverse momentum 
of the parton inside a proton of < kx >= 0.8 — 1.2 GeV. 

Since the Drell-Yan process is expected to be one of the leading back- 
ground contributions at the future high energy facilities, it is important to 
predict its cross section as precisely as possible. The obtained triple dif- 
ferential cross section of the dilepton production in pp collisions is also a 
necessary input for models, studying nuclear medium via high energy dilep- 
tons, produced in pA and AA collisions. To meet this demand and to consis- 
tently evaluate the ISI effects in high energy processes, we need to improve 
our knowledge of the parton spectral function in the nucleon. The single- 
parameter Breit-Wigner parametrization might be insufficient. To pin down 
the parton spectral function, the study of other exclusive processes will be 
necessary. In particular, it should be possible to reduce the sizable uncer- 
tainty in the width. We shall address this issue in the future. 
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A Appendix 

In this appendix we present the details of the Drell-Yan process description 
in different frames of reference. The arguments for the definitions of the 
partonic momentum fractions used in our calculations are given as well. 
Let us consider the Drell-Yan scaling limit (S — > oo). The light cone 
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components of hadron momenta in the center of mass system are 

P?Y = (P?Y = I - Ml ± \l( S -) 2 -M%S. (52) 



Thus, the plus-component of the projectile's momentum P 2 + and the minus- 
component of the target's momentum Pf go to infinity ~ y/S, while all the 
other components are negligible in the scaling limit. 

With the chosen definitions of we get as a limit of (J33|) : 

M 2 = Z&S; 



Applying approximate definition xp ~ 2p z /yS, we recover the well known 
parton model relations: 

zf = 6-6- (54) 

This means that we can use 6 = pf/Pi, £,2= P 2 / P% as the arguments of the 
parton distribution functions in the factorization formula for the Drell-Yan 
process in the center of mass system. 

In contrast, hadron light cone momenta scale differently in the target rest 
frame: 



R± = -±L -M N ± X ( -4r-) - S] Pf = M N . (55) 
2 2M N \ \2M N J 1 v ; 

The motion of the projectile is again confined to the light cone. However, 
there is no special direction for the target parton. Therefore, one needs to pa- 
rameterize the " soft" properties of the target with more general distribution 
functions, which depend on the 4-momentum of the parton instead of a single 
scalar variable 6 These functions W(p) (partonic Wigner distributions) were 
introduced by X. Ji in |H8j . 

M. Sawicki and J. P. Vary [SU] considered the Drell-Yan process in the tar- 
get rest frame within the factorization framework. They used the analogous 
to DIS definitions of the both momentum fractions, 

(Lx=pt/Pi; & = P i/P 2 + ; (56) 
and found scaling violation. Indeed, in the target rest frame (IHHjl transforms 
to 

n/r 2 2 , 2 , £lMN I 2 , S \ , 62P2 / 2 , -2 \ q ~> 

M = mi +m 2 + [m 2 + p 2± ) + —-rr- [rri! + p 1± ) - 2pi±_ ■ p 2± ; 

£2-^2 v > ^iMr 



N 



10 



iiM N &p+ j 
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where we used the definitions (J5fij) and 



Bearing in mind equations and (|57|1. we arrive at the following limiting 
values for M 2 and x^: 




x F = 2£ 2 . (59) 

Hence, the variables (jHfij) do not coincide with the Bjorken variables in the 
Drell-Yan scaling limit. 

Contrary to the statement of jSH], the relations (|59~j) are not a dynamics 
effect, but an artifact, caused by the use of the alternative definitions t)56j). 
At the same time, the factorization in the form applying the usual 
dependent parton distribution functions, is not applicable in this system of 
reference, because the motion of the target parton is not confined to a light 
cone even in the high S limit. 
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